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Abstract 



The dynamics of light fermions propagating in a spatial direction at high temperatures can 
be described effectively by a two-dimensional Schrodinger equation with heavy effective 
mass meff = ttT. Starting from QED, we discuss the transition from three- to two- 
dimensional positronium discussing the latter in detail including relativistic effects. In 
the case of QCD the problem is similar to that of heavy quarkonium. Our effective 
potential contains the usual Coulomb and confining parts as well as a perturbative spin- 
spin interaction. The resulting qq "wave functions" reproduce recent lattice data for the p 
and TT channels. The physical meaning of such 'confinement' is related to the non-trivial 
magnetic interaction of color currents in the quark-gluon plasma. Our results do not 
contradict the idea that the normal electric interaction of color charges is screened and 
produces no bound states in the usual sense. 
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1. Introduction 



It is generally believed tliat hadronic matter at high temperature undergoes a transition to 
a new phase, the so called quark-gluon plasma (QGP), which is a gas of weakly interacting 
quarks and gluons. As far as global properties {e.g., thermodynamic properties) are 
concerned, this statement is supported by both perturbative QCD analyses and numerical 
simulations based on Lattice Gauge Theory (LGT). (See, e.g., general reviews IQ, @).) 
However, this does not imply the absence of non-trivial phenomena even at very high T. 
As an example, one may mention the problem of magnetic screening, the mechanism for 
which remains unknown. 

Some time ago, DeTar [Q, || suggested that the QGP may be 'dynamically confined', 
in the sense that only color-singlet modes produce poles and branch points in linear 
response functions. In other words, he proposed that color-singlet modes might control 
the large-distance behavior of plasma disturbances. DeTar and Kogut P, H pursued this 
matter numerically by calculating the correlation functions of various operators using 
LGT. The 'technical' point of importance for the interpretation of these results is that 
the correlations were measured in the spatial direction. 

These functions were found to decay exponentially with distance, ~ exp(— Mz), and 
DeTar and Kogut ^ studied the temperature dependence of the exponent, M, (known 
as the screening mass) in a number of channels. The screening masses of the chiral partners 
(tt, ct), (p, Ai) and (A^(|+), A^(|— )) were found to become equal within the accuracy of 
the calculations at the temperature where chiral symmetry is restored (defined as the 
point where < qq >^ and referred to below as Tc = T^sr)- 

For spin zero channels, (7r,(T), and in the chiral limit, one would expect the mass to 
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remain small even above Tc. These modes can be viewed as long-range fluctuations in 
the order parameter of the second-order chiral restoration transition as first proposed 
by Hatsuda and Kunihiro p, ^ [1^. Indeed, lattice data show that M^(= M^-) remains 
relatively small until rather large T. (Further discussion of this point and a compilation 
of data can be found in [|ll|].) In contrast to this result, Eletskii and loffe [jl2[ noted that 



for the other doublets the following is true: Mp = M^^ — 27rT = 2ujq. For the baryonic 
doublet the mass is about SuT = Sujq. These values are essentially those suggested by 
the lowest quark Matsubara frequency, = ttT, and led to the suggestion that these 
correlators can be interpreted as describing the independent motion of two (or three) 
quarks in a plasma described as an ideal quark gas. 

For people not comfortable with the Matsubara formalism, one can explain this point 
in a different way by looking at equal time correlations, as recently suggested by one of 
us [|13|. One can show that the propagator for free massless quarks at temperature, T, 
over a spatial distance x in zero time is given by 



St{x) = (7^9^) J ,3^ exp(ifc-f)[- - - 



(27r)3fc ''2 l + exp(fc/T)' 

The last term contains Fermi occupation numbers, and the combination (1/2 — npik)) is 
clearly related to the familiar combination {1/2 + nsik)) in the case of Bose statistics. 
In both cases the '1/2' term is due to zero-point oscillations and the last term to the 
contribution of thermal excitations. However, quarks with small momentum have Uf = 
1/2. Thus, one has a cancelation of these two contributions for A; <^ T which eliminates 
power contributions and leads to St{z) ~ exp(— vrT^). Another amusing way to look at 
this comes by rewriting (1/2 — np{k)) as (1/2)[(1 — up) — np] so that the particle-hole 
symmetry np ^ {1 — rip) is made obvious. Since the situations near n/ = and n/ = 1 



are physically similar (up to the sign of the propagator), it is natural to expect that the 
midpoint, np = 1/2, is special. There, the particle and hole terms cancel. 

It is useful to mention that there are essentially two simple ways to make the connec- 
tion between the imaginary-time formalism and real time. One is to look at equal time 
correlations {i.e., correlators integrated over frequency) as in the previous illustration. 
The other possibility is to study correlators at zero frequency {i.e., correlators integrated 
over imaginary time.) It is the latter which is used to extract the screening mass in LGT 
calculations while the former is used to determine the wave functions. 

The scenario of free motion in a ideal quark gas is too simplistic to describe the real 
situation, but the observation that most screening masses are consistent with this picture 
suggests that other mechanisms have little net effect. In particular, these data provide an 



important limitation on any dynamically generated quark mass, nidyn [|14[. The screening 



mass for massive quarks would he M = ^ uo^ + rn^y^. Thus, it seems that an appreciable 
rridyn is excluded even in the region T ~ T^sr- A picture in which hadrons go massless 
1^ as T — > Tj^sR (leading to the scenario suggested in [l^, Q) is not excluded. (However, 



given a temperature of T ~ 150 MeV, an accuracy of about 5% for the calculated screening 
masses is required to rule out dynamical quark masses smaller than rridyn = 150 MeV.) 

The next important step was made by recent lattice studies ||T^ which revealed that the 
qq pair, moving over large distances in the spatial direction, is actually closely correlated 
in the transverse plane. Moreover, this correlation is even more pronounced at finite 
T ^ 150 MeV than it is at T = 0. It was suggested \Tf\ that, in line with DeTar's 



dynamical confinement ideas, the quark and antiquark are not free at all but bound even 
at high T. The aim of the present paper is to provide a quantitative explanation of these 
data as well as a general discussion of what these data actually mean for the physics of 



the QGP. 

An essential technical point, important for the understanding of the terminology used, 
is the interchange of the time and z-axis which we shall describe in some detail in the 
next section. After transformation to this 'funny space', one deals with a system at 
zero temperature but placed in a box periodic in the 2;-direction with periodicity f3 = 
1/T. Using this language it is obvious that, in the high temperature limit, a so-called 
dimensional reduction takes place. The (3 + l)-dimensional gauge theory becomes a 
(2 + l)-dimensional theory as suggested long ago |T^. The main physical point is that, in 
the high T limit, the motion of a quark in the 'funny space' is dominated by its momentum 
in the 2;-direction which is governed by ttT as a consequence of the antiperiodic boundary 
conditions. For motion in the transverse direction, as we shall show, this momentum 
behaves like a mass. Since at high temperature this 'mass' becomes very large, any 
attractive potential can bind quarks in the transverse direction. 

The main objective of this work is to derive the high-T limit of the "wave functions" 



introduced (and measured on the lattice) in |]T^. In section 2 we start with the QED case 



and discuss how the transition from d = 3 to d = 2 positronium takes place. Qualitative 
features of = 2 bound states in the d = 2 logarithmic Coulomb potential were recently 
discussed in ref. [113. We report the results of a more quantitative analysis including next- 



order corrections. We pay special attention to the spin-spin interaction because this issue 
is important for the mesonic splittings discussed in section 3. 

The effective d = 2 QCD potential for a gg-pair propagating in the spatial direction 
was studied on the lattice in ref. [Q. It was found to be essentially temperature inde- 
pendent. (This result stands in sharp contrast to the usual potential, for propagation 
in time, which exhibits both deconfinement and charge screening.) Although the data 



pO| are not very good at large separations and allow a wide range of interpretations, 
comparison with the "wave functions" allows us to conclude that the "string tension" 
is essentially T-independent in the interval between T = and T = 210 MeV studied 
in ref-ljl^. Another issue addressed in section 3 deals with the differences, both in the 
shape of the "wave functions" and the screening masses, between the vr and p systems. A 
spin-spin interaction of the Fermi-Breit form provides a natural explanation of most of 
these differences. 

In section 4 we proceed to a discussion of the physical meaning of these "bound states" . 
In order to understand these results, one must come back from the 'funny space' to ordi- 
nary space. One realizes that by dealing with static correlation functions (with no time 
involved), it is not possible to say anything about the energy spectrum of excitations in 
the QGP at high temperatures. What is actually studied is the momentum spectrum of 
the correlators for various channels. At T = the two spectra are identical. At high T, 
they are completely different. 
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2. Propagating an e+e pair at high T 



In this section we shall study the "wave function" of an e^e~ pair at finite temperature 
in Euclidian space. We shall proceed by ignoring the Coulomb force and restricting our 
attention to the magnetic interaction between currents. The aim of this restriction is to 
preserve the analogy to QCD in which the electric interaction is screened and magnetic 
interactions are believed to be dominant at high temperature. The purpose of this example 
is that QED is an abelian gauge theory with a fixed and small coupling constant, e. One 
can work out all essential ingredients, and both the "binding energy" and the "wave 
function" of an e+e" pair at high T can be obtained. As is well known, a quantum field 
theory at finite temperature, T, can be formulated using Euclidean time, r, ranging from 
to 13 = 1/T with periodic boundary conditions for bosons and antiperiodic boundary 
conditions for fermions. This formulation, which is the corner stone of LGT calculations, 
will be the starting point for our considerations. This convenient generalization of the 
T — case has, however, its limitations. One is unable to consider any time-dependent 
quantities and deals only with static (heat bath averaged) quantities. As a consequence, 
LGT cannot tell us anything about the frequency spectrum of elementary excitations. 
What we are going to study in this work is, therefore, related to the momentum spectrum. 
(In other words, properties related to spatial correlation functions.) 

Given the magnetic nature of this problem, it is simpler to proceed by interchanging 
the (compressed) Euclidean time, r, and the spatial direction, z, in which the correlator 
is measured. This transformation leads us to a 'funny space' where the system propagates 
in unlimited time (hence, at zero temperature) but is confined in a box in the 2;-direction. 
This transformation must be applied to the fields. Specifically, the original magnetic field, 
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Fi3 and F23, becomes an electric field {i.e., F13 =^ Fio = and F23 =^ F20 = Ey). The 
magnetic interaction between two rapidly moving charges is transformed into an electric 
interaction between two, nearly static, charges. We emphasize that this new Coulomb 
interaction has emerged from a technical manipulation. Its physical origin remains mag- 
netic and, as such, it is not subject to the effects of screening which are presumed to 
suppress genuine electric interactions. 

The presence of a spatially periodic "box" modifies the interaction of the charges. 
This modification is readily calculated for our simple example of small (a = <^ 1) 
charges in QED with dynamics given by the current-current interaction of the leading one- 
photon exchange diagram. The r;z-interchange immediately transforms this interaction 
into the ordinary Coulomb potential, V — —a/R. The presence of the "box" boundaries 
modifies this result significantly because the electric potential must now satisfy periodic 
boundary conditions in the ^-direction. There are two ways of solving the standard 
Poisson equation: 

- V^* = p. (1) 

One can either (i) expand the field in Fourier components 

n 

p(a:,y,;s) = ^p„(a:,y)e-'2n.T. ^2) 

n 

(in the high T limit only the n=0 modes are important) or (ii) one can construct a periodic 
array of "reflected charges". (See Figure 1.) The second approach is certainly simpler; it 
involves no integrals and leads directly to the sum 



We note that the sum is automatically regularized at large n and that the large-T limit 
of this three-dimensional problem is simply the two-dimensional Coulomb potential: 

which is equivalent to the potential of a uniformly charged wire (with charge density e//3) 
in three dimensions. 

The fermion field is also affected by the T2;-interchange. This is most easily appreciated 
by considering the Dirac equation for a single, massless particle in Euclidean space which 
has the form: 

l^d.iP = (5) 
where the 7^ are Euclidean gamma matrices. The solutions to this equation are 

ip{x) = ^06"^^^"+^^''^''° e-'P-^ (6) 

where the dependence on the 0-direction is due to the antiperiodic boundary conditions 
which result from the finite temperature. (Here, iJjq is a spinor which will not concern 
us.) The Euclidean Dirac equation for this problem is evidently invariant under the tz- 
interchange, and all effects are concentrated on the wave function, eqn.(6), which becomes 

because of the boundary conditions. The Dirac equation is solved by the choice 



po = ±i VSTpi = ±i ^J {{2n + l)nTf + pI = -lE (8) 
with E real. The wave function thus takes the form 



which is the solution of a Dirac equation for a particle with a given momentum, pz — 
(2n + l)7rT, in the z-direction. For large Xq only the state with lowest energy, E, survives. 
This state corresponds to the ground state of the associated real time Dirac equation 



7^i?^ = (7f P± + 7f (2n + l)7rT) ^ (10) 

where the superscript, M, now indicates Minkowski 7-matrices. This equation is equivalent 
to a Dirac equation in (2 + l)-dimensions with a "chiral" mass of m = (2n + l)7rT. This 
can be rotated into a real mass by the unitary transformation 

Thus, we have 

-i^Ei; = (7f + {2n + l)7rr)V'. (12) 

So far we have discussed the case of a massless particle. For massive fermions, one should 
use an effective mass of 



^eff = + [(2n + l)7rTp. (13) 

Since lattice gauge calculations provide no conclusive evidence of a dynamical mass, we 
will restrict ourselves to the case of vanishing dynamical mass. 

It is clear from eqn.(13) that the effective mass of the electron is very large at high 
temperature. The dynamics thus become nonrelativistic and can be described by a 
Schrodinger equation in two dimensions: 

2 

Ejjj^^ijj. (14) 
^ 2m ^ ^ ^ 
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The effects of any additional potential {e.g., an "electrostatic" potential, V = — e$) can 
then simply be included by making the substitution E ^ E — V in eqn.(]T^). 

These manipulations are readily extended to the present example of an interacting 
e~^e~ pair. (The details are given in the Appendix.) We note that the restriction to zero 
Matsubara frequency, n = 0, in eqn.(|35|) represents a decision and not an approximation. 
It is also desirable to make the approximation of neglecting all terms, in eqn.(|35D with 
k ^ 0. As we shall indicate shortly, this approximation introduces a tolerable relative 
error of O(e^) for large T. 

Before turning to numerical results, it is useful to perform a dimensional analysis 
aimed at understanding the high T behavior of this system. For a wave function of size, 
R, the kinetic energy is equal to (vrTi?^)"^. Neglecting the binding energy, E, we obtain 

-^-^ ^ e^T. (15) 

Thus, the size of the bound system and the binding energy are 

~ — and E ~ e^T. (16) 
eT 

The size of the system is identical to electric screening length in a hot plasma which is 
also (eT)-i . § 

With the approximations indicated above, the reduced mass for the e~^e~ pair becomes 
m = vrT/2, and the Schrodinger equation, eqn.(^), assumes the form 

Using dimensionless units, 

pr 

(18) 
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we obtain for the / = state 



( 



(19) 



where 



-kE 



+ 27r ln(2Ta) 



(20) 



E 



eT 



and 



Veffix) = 271 ln(x). 



(21) 



Note that eqn.(|T9D no longer depends on the physical parameters. Here, we have used the 
high temperature limit of the effective potential given by eqn.(^). The wave function is 
now given in terms of the dimensionless variable, x, rather than the physical parameters, 
e and T. We have solved eqn.(|l^) numerically and the resulting wave function is plotted 
together with the potential in figure ^ Restoring dimensions, the energy is 



We now consider a variety of corrections to these high-T results which, in the approx- 
imation of zero dynamical mass, can be expressed purely in powers of our dimensionless 
parameter, the electric charge e. In particular, we consider (i) the non-asymptotic part of 
the potential, (ii) excitation of fermion modes with higher relative Matsubara frequencies 
(/c 7^ 1), and (iii) relativistic corrections. 

In figure 3 we compare the modified Coulomb potential in the periodic box given by 
eqn.(§) with its high-T, two-dimensional limit for a few values of the temperature. One 
can see that, although agreement is generally good enough, there are significant deviations 



E = e2T(1.45 -21n(e)). 



(22) 
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at small distances: The original Coulomb potential, e'^/Anr, is still present in the original 
form. Since a typical separation is r ~ 1/eT, this produces corrections to the energy of 
the system of O(e^T). 

The approximation of neglecting all terms in eqn. (|35|) with k ^ can be relaxed 
perturbatively. Applying standard second-order perturbation theory one obtains 

^E=-T.^^}^^-e^T. (23) 

Estimating the matrix elements as V ~ 0{e^T) and all energy denominators as 0{T), 
one finds this correction to be 5E ~ e^T. As expected, this correction is smaller than the 
leading contribution to the energy, eqn. (25), by a factor of O(e^). Taking only the first 
excited {k = 1) state into account numerical integration has produced the following result 
for = 0.2 : 

SE , 

— = 3.5 X 10~^ (24) 
E 

A typical velocity, Vt, transverse to the compact 2;-axis is Pt/f^^es ~ l/RT ~ e. (The 
somewhat surprising feature of this result is that it is not as it would be for the d = 3 
hydrogen atom.) Therefore, relativistic corrections start at the level i^rei ~ ^^f^-Enon-rci ~ 
e^T. A relativistic expansion similar to that leading to the Fermi-Breit Hamiltonian can 
be performed as in the d = 3 case. We again have —p'^/Sm'^ kinetic energy corrections, 
the spin-orbit term (which does not affect the ground state), spin-spin and "Darwin" 
terms. 

We shall be concerned only with the properties of the spin-spin interaction since we 
wish to discuss the spin splitting in QCD in some detail in the next section. In the d = 2 
space one has symmetry only with respect to rotation in the xy-plane. Hence, states 
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should be characterized by the z-component of their spin. In the ground state (with 
vanishing 2;-component of orbital angular momentum), the e^e~ pair has two degenerate 
states with total spin 5*^ = ±1 and two distinct states with 5*^ = of different symmetry 
representing total spin S = 1 and 5 = (ortho and para-positronium, respectively). 

We start with with the textbook {d = 3) expression for the Fermi-Breit interaction 
between two opposite charges 

Vsreu = ^[-^ + + y^i ■ <^AR)] (25) 

which consists of the familiar "tensor force" between magnetic moments and a contact 
interaction. In addition, in the case of positronium, one has to include an "annihilation" 
contribution^ 

Vann = :^{S + a,-a2)6{R). (26) 

Starting from these formulae, we can obtain the d = 2 effective interaction (needed in 
the high T limit) by averaging these potentials over z-direction. In three dimensions 
the tensor force vanishes for the spherically symmetric ground state. However, for finite 
temperature, the ground state has a "pancake" shape and a non-zero quadrupole moment. 



This leads to a contribution from the tensor force. This can be seen by rewriting eqn. (^5D 
in momentum space as 

VBreit = -. 5- / xn e ^ [0-1-0-2 5 J- (27) 



Am? J (27r)3 g2 

Averaging over z means that we select the component of the potential which has no 

momentum transfer in the z-direction. Therefore, the second term contributes only if the 

•^In the next section (in which we discuss the quark-antiquark interaction), the annihilation term is 
not present. 
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spins are in the x?/-plane:0 



I / dZ VBreit = o / /o \o ^ ' ^2 



4m2 y (27r)2 ^ J 



^It is useful to consider an alternate derivation of the spin-spin interaction which offers a different 
explanation for why spins in the z-dircction represent a special case. As for the Coulomb field in the box, 
one can satisfy periodic boundary conditions for the magnetic field by using a periodic array of reflected 
magnetic moments and by then noting that the z-directed moments display a strong tendency to cancel. 
A nice way of explaining this was suggested to us by H.A. Bethe. The Coulomb potential is a coherent 
sum of all reflections as they are all of the same sign. At high T, the potential of a charged wire results. 
Magnetic moments, on the other hand, can be considered as a pair of magnetic monopoles sitting on this 
wire with alternating signs. From this it is clear that the magnetic interaction vanishes since, for large 
T, the distance between these monopoles decreases. 
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The remaining local terms simply lead to two-dimensional 5-functions. 

\ f dz8\r)=T5\n). (29) 
p Jo 

As a result, the spin splitting of d = 2 positronium is 

_ WT|^(0)P 

^spin—spin ^2 ^ v / 

with the spin factors C„ = 3/4,1/2,-1/2 for 5' = 1,5*2 = =tl) ^S* = 1,5*2 = and 
5 = 0, = states, respectively. 
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3. qq "binding" in hot QCD 

As noted in the introduction, there exists a significant general difference between QCD and 
QED. Even at arbitrarily high temperature, one is confronted with non-trivial manifesta- 
tions of the non-abelian nature of the former theory as one can see from the "dimensional 



reduction" argument [0. The resulting d = 2 + 1 theory is far from simple. In par- 
ticular, as argued in |^T], |22|, it is known to produce a linear confining potential just as 
d = 1 + 1 and = 3 + 1 (real QCD) theories do. In QED the electron charge is a physical 
parameter, and this can be assumed to be small in order to make the perturbative theory 
meaningful. In QCD we have a "running" coupling constant, and the actual magnitude 
of Us = g'^/4:TT depends on the distances (or momenta) relevant for a given application. 
The precise value can only be specified if higher-loop corrections are calculated. In the 
high-T limit, for example, we have a large effective mass and therefore the size of the 
Coulomb bound state is small, R ~ 0{l/gT)!^ However, the logarithmic corrections are 
actually cut off at the smaller scale of 1/T. Therefore, (in contrast to quarkonium) one 
should expect that the relevant as is as{T) which becomes small as T ^ oo. 

Unfortunately, the lattice data which we wish to discuss do not correspond to this 
relatively clean limit because they have been obtained at temperatures close to the critical 
temperature {i.e., T ~ 1.5Tc ) which corresponds to effective masses of rries = ttT ~ 
450 — 600 MeV. From the experience of J/ ip and T physics, we know that these masses 
are not sufficiently large to justify use of the Coulomb force alone. However, we have also 
learned another lesson from these examples: A non-relativistic approach based on a more 



^This is similar to the case of superheavy quarkonium pSf with the obvious substitution 
(4/3)(5V47r). 
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complicated effective potential (which includes confining forces) works well in describing 
quarkonium bound states. It is reasonable to try the same approach in our "not-so- 
high" T case. In practice this means that, as in the previous section, we will rotate 
into the 'funny space' and solve a Schrodinger equation for quarks of mass rrics = ttT. 
Consequently, the effective potential must also be given in this rotated space. Fortunately, 
such an interaction has been extracted on the lattice at high T for pure glue pO|- The 
results are shown in figure ^. Manosakis and Polonyi fit their results using the standard 
parameterization 

V = —a/r + ar + constant (31) 

with a = 0.184 ± 0.02 and ^/o^ = 0.22 ± 0.03. These results can be compared with 
the values a = 0.26 and y/a = 0.22 ± 0.02 obtained for T = 0. If one identifies a 
with the phenomenological value of the string tension, (400 MeV)^, the physical size of 
the lattice spacing is a ~ 1/(2 GeV) ~ 0.1 fm and T ~ (1/3) GeV. The constant 
in eqn.(^) was not evaluated in ref.[^. It is also instructive to compare these data 



with the potential appropriate for quarks propagating in the time direction which is 
also shown in figure ^. Not only are the confining forces absent in this case (because 
we are well above the deconfinement temperature), but the Coulomb potential is much 
smaller due to screening effects. In contrast to these results, the potential measured in 
the spatial direction appears to be essentially T-independent, apart from some change 
at small distances. (This difference may simply result from the decrease of the effective 
Coulomb charge due to a difference in radiative corrections as g{r) g{f3)-) 

At this point one may ask how much of the LGT potential is simply due to the modified 
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Coulomb potential obtained in the previous section, 



V = a{-l/r + 2T\og{rl[3)) (32) 

which results from the inclusion of "reflections" of the charge in the box. Fixing the 
coupling constant at small r, one obtains the curve plotted in figure H which is not excluded 



by the LGT data. Thus, the data of ref.||20|, do not rule out the possibility that the "spatial 
string tension" is actually very small. However, at the lower temperature, T ~ 210 
MeV, where the 'wave functions' are being measured, the difference between the modified 
Coulomb potential and the string potential are substantial (under the assumption that 
the string tension remains the same). Therefore, by studying the wave function one can 
discriminate between these two alternatives. As we will show in the following, the 'wave 
functions' for vr and p given in ref. |]T^ seem to require the presence of a string potential 
of the above magnitude. 

With the potential fixed, one can proceed further and solve the d = 2 Schrodinger 
equation numerically. In figure ^ we show the corresponding wave functions for sev- 
eral temperatures for both potentials eqn.(|31|) and eqn.(|32|). These results are compared 
to data at T = 210 MeVQ These results suggest that the confining effective potential, 
eqn.(0), reproduces the data well while the modified Coulomb potential is definitely too 
weak at large distance^. 

We now consider the role of spin-dependent forces. Comparing data for the p and vr 

channels (shown in figure ^), one notes significant differences. The average size of the pion 

^The data are taken at /3 = — 5.445 with four time slices. Using the same (3 and quark mass but 
six time shces, the authors derive a temperature of Tg = 140 MeV This leads to a temperature of 
T4 = 210 MeV for four time slices. 

^Provided the coupling constants for the logarithmic and 1/r parts of the potential are the same. 
Indeed, in this case, the standard 'asymptotic freedom' arguments do not apply. Logarithmic corrections 
are cut off at a scale of 1/T, so both terms should have the same coefficient, as{T). 
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is significantly smaller than that of the p. At small distances (r < 2 a), the pion wave 
function is much largeiQ with |'?/'(0)p being about four times greater than in the p-meson 
case. 

It is natural to ask whether this relative compactness of the pionic gg-wave func- 
tion results chiefly from the perturbative spin-spin interaction due to the exchange of 
high-frequency magnetic fluctuations. The corresponding (z-averaged) interaction was 
discussed at the end of the previous section, and it can simply be used after (i) switching 
from standard QED to QCD units, (4/3)((7^/47r), and (ii) ignoring the annihi- 

lation contribution in the case of these isospin 1 channels. Because the effects of this 
spin-spin interaction are not small, we cannot simply estimate them using perturbation 
theory. Rather, we must re-solve our effective Schrodinger equation including this spin- 
dependent, local interaction]^. 

Two sets of results at temperatures of 210 and 350 MeV are shown in figure ^(a) and 
(b). The three curves correspond to {S = 1, = ±1), {S = 1, 5"^ = 0) and {S = 0, = 0) 
states. A value of the effective coupling constant of = (7^/471 = 0.20 was required 
to describe the difference between tt and p wave functions found in LGT calculations. 
This value of is about 50% larger than that used in the central interaction, eqn.(^Tp. 

This difference may have several origins, but should not taken too seriously because, at 

^It is worth noting that the T = pion wave function measured by Bernard et al. is in principle 
the same as the zero temperature Bethe-Salpeter amphtude calculated by Chu, Lissia and Negele psf . 
Both drop with an initial e-folding scale of 0.48 fm which would give the pion an rms-radius of ~ 1/3 
fm. This is substantially smaller than the charge radius obtained in lattice gauge calculations of the 
density-density correlation function . In terms of phenomenology ||2^ , the Bethe-Salpeter amplitude 
provides the "intrinsic size" of the pion (probing its gg-component) while most of the charge radius 
actually comes from the coupling of a virtual photon to the p-meson "cloud" . 

^In order to avoid difficulties arising from a zero-range potential, we represent the (5-function appearing 
in VBreit by 5{x) = sin(A2:)/7rx. Here, A represents a high frequency cutoff which we have chosen to be 
A = 7r/a in order to simulate the effects of a finite lattice spacing, a. 
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temperatures of T ~ 210 MeV, the wave functions are rather insensitive to the strength 
of the central Coulomb potentiaQ We remark that the lattice wave functions are not 
normalized. In figure |^(a) the two different states of the p have the correct relative 
normalization, while the overall normalization of vr and p-wave functions is arbitrary. In 
figure ^(b) all three wave functions have the same normalization. 

Our final point concerns the splitting of the screening masses. From the calculation 
displayed in figure we obtain 

Ep-E^ = 340 MeV (33) 
which can be compared with the results of direct measurements on the lattice with six 



time slices {i.e., at the lower temperature (T = 140 MeV), which gives ~ 440 MeV. [p^ 

These predictions can be tested much more accurately by looking at the difference 
between the longitudinally and transversely [Sz = ±1,0) polarized p or by going to 
higher temperatures where our d = 2 approximations should be better justified. Although 
the differences between wave functions for the different states of the p shown in figure 
^ appears to be rather small for both temperatures, their ratio (particularly at short 
distances) should provide sufficient information and should be less affected by systematic 
errors. The screening mass of these states differ by AM = 0.25T 

The general question of the physical origin of spin-dependent splittings in hadronic 
physics has been much debated in the literature. Together with the Breit-type perturba- 



tive interactions, interactions due to instantons, as proposed in refs. |3y, ^ , also produce 
a quasi-local interaction of the form (c?i ■ (?2)(ii ■ ^2)- O'^^ distinguish them is to 

consider their temperature dependence, because the instanton density is expected to de- 

i°We point out that at T = Karsch |2| gives a value of as{Tc) = 0.29 for pure glue. Also, the T = 
value of as = 0.3 extracted from charmonium decay is of the same magnitude. p9| . 

21 



crease rapidly at high T, while the perturbative mechanism can even be enhanced as the 
bound states discussed above become more compact. 

Finally we mention that the effect of the Coulomb piece of the interaction, eqn. (|3ll) , on 
the wave function is tiny at those temperatures where the lattice data have been obtained. 
Therefore, it seems that these wave functions require the presence of a string-like confining 
potential. 
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4. Conclusion and Discussion 



Before summarizing our results, let us discuss the physical meaning of the spatial cor- 
relation functions. It has been already mentioned instead of yielding the more familiar 
spectrum of excitation energies, spatial correlators provide only the spectrum of momenta 
(or screening masses). In distinction to the situation at zero temperature, where these 
spectra are identical, the energy and momentum spectra are by no means similar at 
finite temperatur^^ Unlike the former, momentum spectra display a very smooth T- 
dependence. There are two reasons for this fact. First, the quark effective mass consists 



of two parts, rries = y vr^T^ + m'^y^. The first increases with T and the second decreases 
as T ^ Tc with the result that the total value is roughly constant below Tc. Second, due 
to some dynamical reason (as yet unclear), the effective potential for spatial propagation 
is also temperature independent^. As a result, there is a correspondence between the 
poles (and the corresponding "wave functions") for the T = and high-T cases. In this 
(but only in this) sense, hadronic degrees of freedom smoothly traverse this boundary. 

It is also worth recalling that the trend seen in the momentum spectra discussed above 
is quite different from what is known about the trends of the energy spectra. Various 
theoretical models suggest that the energy gaps (usually still referred to as "hadronic 
masses") decrease with T. (See, e.g., recent calculations based on hadronic scattering 

data [^.) There are suggestions that they even vanish at T = |T3| so that, at higher 

^^We mention one more simple example which emphasizes the difference between the two. Consider 
non-interacting, massless fermions at high T. The excitation energy spectrum obviously starts from zero 
and corresponds to the cut while the momentum spectrum studied is a sequence of poles at (2n + 1)t:T. 
The lowest pole is our "quark effective mass" in the 'funny space'. 

""^^Apparently, this potential "does not notice" even a very strong first-order deconfinement transition 
which is observed in pure giuonic theory. Notice that this potential has not yet been evaluated when 
dynamical quarks are present. 
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temperatures, the energy spectrum is just a cut containing excitations of any energies 
starting from zero. There is no contradiction between these two pictures: Propagation 
in space and in time at high T are no longer simply related. Each can tell us something 
interesting about the properties of the QGP. 

In this paper we have addressed not only these momentum (or "screening mass" ) spec- 
tra but also the corresponding wave functions. We have shown that the spatial correlation 
functions at high temperatures can be quite well understood in a 'funny space' obtained 
by interchanging time and space axes. The two main ingredients of our calculations are 
the quark effective mass, frieff = ttT, and the effective potential, Veff, which is a combi- 
nation of a Coulomb and a confining part and is supplemented by a short-range spin-spin 
interaction. We find these ingredients to provide a good description of the wave functions 



measured in |T^. We further find that, at the temperatures for which lattice wave func- 
tions have been measured (T ~ 210MeV), the effect of the Coulomb piece of the effective 
interaction is small. Hence, we conclude that the lattice data require a relatively strong 
potential — much stronger than the modified Coulomb potential which would arise from 
dimensional reduction arguments [|l9l. Thus, as expected, dimensional reduction is not 



valid at temperatures as low as T ~ 210 MeV. Additional lattice calculations at consid- 
erably higher temperatures would be needed in order to check the validity of dimensional 
reduction. 

The issue certainly can be pursued further both numerically and analytically. In 
particular, better data and data at higher temperatures are needed for further tests of 
the theory, as well as more detailed studies of the effective potential by itself a la pO| . 
For example, calculations of the effective potential with dynamical quarks are needed. 
Obviously, further attempts can be made in order to get a smooth matching between the 
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traditional potential model description of hadrons at T = and our approach at high 
T. This could be done by solving the complete three dimensional Dirac equation with 
anti-periodic boundary conditions. Splittings between different spin-isospin components 
can be studied in greater detail. 

What have we actually learned from these "wave functions" in simple physical terms? 
Consider a spatial correlation function at distance, z, with the two ends taken for sim- 
plicity at the same instant of time, and let us describe qualitatively the typical paths 
contributing to it. (This language is convenient because it is the same in physical and 
'funny' space.) At T = the light quarks may travel by any pathQ deviating from the z 
axis hj 6t ^ 6x ^ 6y ^ z (see figure 0(a)) which leads to power-like correlation functions. 
At high T, for non-interacting quarks, the deviations in time are smaller, 6t ~ l/T, due 
to the cancellation effects described in the introduction, and the correlator is exp(— Mz). 
As a result, 6x ^ 6y ^ \J z/T. (See figure 0(b).) These are the paths of "heavy" non- 
interacting quarks in the 'funny space'. Finally, for high T and interacting quarks, quark 
and antiquark now travel along paths which are correlated within a distance governed 
by the "wave functions" considered above. (See figure 0(c).) Their deviation from the 
straight line remains as in case (b). 

The perturbative Coulomb interaction in 'funny space' arises physically from the per- 
turbative magnetic interaction of two color currents in original physical space. In the very 
high T limit, when this effect dominates, the distance between them is (5x ~ 5 ~ l/{gT). 
Note that in a quark-gluon plasma the magnetic interaction is screened (if at all) only on 

a parametrically larger scale of Rm ~ ^/{q^T). Thus, there is no contradiction here. 

^■^ As usual, one can say that the correlation is caused by pair production somewhere in the middle, 
and its existence does not contradict causality because one cannot use it for signal transfer. 
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Of course, the origin of the non-perturbative "confining" potential remains unclear. 
It is, in any event, due to some non-perturbative Gxy field which is the magnetic field 
both in physical and 'funny' spaces. 

It is important to distinguish between magnetic and electric interactions which are 
very different even in the more familiar electromagnetic plasma. The electric fields are 
screened and appear only as high frequency fluctuations. Magnetic fields are not, and 
they can produce very spectacular structures (as in the sun). In QCD the additional 
non-trivial feature is that the long-range magnetic fields are self-interacting, so their 
structure at high T is not understood at all. 

As a final remark, let us try to answer the following question: Are the phenomena 
considered here subject to experimental observation? Correlations of electromagnetic 
currents are the source of photon and dilepton production from the plasma. Therefore, 
correlators similar to those introduced here will be essential in determining the effects of 
perturbative and non-perturbative gluomagnetic fields on these reactions. Such correla- 
tors are not easily obtained from LTG calculations. Here we have shown that a simple 
dynamical model is capable of reproducing the results of these more complicated LGT cal- 
culations for zero frequency correlators. This calculations have established a consistency 
between the zero frequency correlator and the spatial potential of ref.[^. Moreover, they 
reveal the absolute necessity of a linearly rising piece in this potential. Given the success 
of this dynamical model and its relative ease of implementation, it seems both practical 
and promising to use it for more detailed calculations of photon and dilepton production. 
Thus, although much additional work needs to be done, it may turn out that the topic of 
this paper is not as academic as it may appear at first sight. 
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Appendix 



In this appendix we will describe the details of the Schrodinger equation for an e~^e~ 
pair. The wave function is characterized by the Matsubara frequency of the pair in 'funny 
space', Pz = 2miT. The frequencies of the electron and positron must sum to this value 
but their difference is free. The wave function for the e^e~ pair can thus written as 

oo 

^'^= E c.^fc (34) 

fc=— oo 

where n is the sum of the Matsubara frequencies and — oo < k < +00 their difference. The 
coefficients, c^, are determined from the solution of the following Schrodinger equation: 

(-^'l" + '= + ^l + l"-*-^l' + 2(.n+f+l).r + 

„2 \ +00 

_^l^_^V„y^^^J^V.Ai.t..UU = E^i.t (35) 
where Vo is the effective Coulomb Potential (eqn.(|^)) and 

are the higher Matsubara components of eqn.(|[). 

These energies and wave functions are the immediate building blocks for the various 
correlators which we can construct. In the original, unrotated space, these are given by 

C{f±,z,T) = J d'^x± < T[0{x±, z = 0,r = 0)q{x±, z,T)Tq{x± + r±, z,t)] > 

= '^exp{-iuJnT)C{f±,z,Un) (37) 

n 

where the "source" O is an operator which creates quark and antiquark at z,t = 0, F 
the appropriate Dirac matrix for the meson under consideration and Un are the usual 
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Matsubara frequencies. The Matsubara components of the above defined correlation 
functions are then related to the solution of the Schrodinger equation (|35| ) by 

C{r^,z,uJn) = 5:a,exp(-M^z)^^(r-l) (38) 

a 

The index, a, labels the eigenstates of eqn.(p5D, are the associated screening masses 
of these states which are identical with the eigenvalues of the Schrodinger equation (|35|) 
and Qa are some coefficients which we do not further specify. In the limit of large z only 
the ground state survives and we have 

Jim C(r^, z, Un) cx exp(-Mo"2)^^(r^) (39) 

where '^q denotes the solution of eqn.(^) of lowest energy for given z-momentum in 
'funny space' Pz = ^unT = ojn- 

The primary focus of LGT so far has been the zero frequency correlator {i.e., n = 
0). This is achieved by integrating the correlator, eqn. (p7|) , over r. In addition, an 
uncorrelated source operator of the form 

Ouncorr.ix±, Z,t) = J dr' (fx'^ q{x±, Z, t) T q{xi_ + Xj_, Z,T + t') (40) 

is used. This source sends off a quark-antiquark pair with vanishing relative Matsubara 
frequency and transverse momentum. 

Finally, notice that transverse wave functions result from correlations between quark 
and antiquark at equal time because both operators are defined at the same time, r, in 
eqn.(|3). 



29 



References 

[1] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 

[2] E. Shuryak, Phys. Rep. 61, 72 (1980). 

[3] C. DeTar, Phys. Rev. D 32, 276 (1985). 

[4] C. DeTar, Phys. Rev. D 37, 2328 (1987). 

[5] C. DeTar and J. Kogut, Phys. Rev. Lett. 59, 399 (1987). 

[6] C. DeTar and J. Kogut, Phys. Rev. D 36, 2828 (1987). 

[7] K. Born, S. Gupta, A. Irback, F. Karsch, E. Laermann, B. Petersson, and H. Satz 
(MTe Collaboration), Phys. Rev. Lett. 67, 302 (1991). 

[8] T. Hatsuda and T. Kunihiro, Phys. Lett. B 145, 7 (1984). 

[9] T. Hatsuda and T. Kunihiro, Phys. Rev. Lett. 55, 158 (1985). 
[10] T. Hatsuda and T. Kunihiro, Prog. Theor. Phys. 74, 765 (1985). 
[11] A. Gocksch, Chiral symmetry in hot QCD, Brookhaven preprint BNL-46286,1991. 
[12] V. Eletskii and B. loffe, Sov. J. Nucl. Phys 48, 602 (1988). 

[13] E. Shuryak, Invited talk at Quark Matter-91, Gathnburg, 1991, Nucl.Phys.A, in 
press. 

[14] H. Bethe, G. Brown, A. Jackson, and P. Pizzochero, subm. to Nucl. Phys. A. 

[15] G. Brown and M. Rho, Phys.Rev.Lett. 66, 2720 (1991). 

30 



[16] G. Brown, H. Bethe, and P. Pizzochero, Phys. Lett. B 263, 337 (1991). 

[17] C. Bernard, M.C Ogilvie, T. DeGrand, C. DeTar, S. Gottlieb, A. Krasnitz, R. Sugar, 
and D. Toussaint, 'The spatial structure of screening propagators in hot QCD' 
University of Arizona report AZPH-TH/91-60, (1991) 

[18] T. Appelquist and J. Carazzone, Phys. Rev. D 11, 2856 (1975). 

[19] T.H. Hanson and I. Zahed, Stony Brook preprint SUNY-NTG-91-44. 

[20] E. Manosakis and J. Polonyi, Phys. Rev. Lett. 58, 847 (1987). 

[21] E. D'Hoker, Nucl.Phys.[FS4] B 200, 517 (1982). 

[22] C. Borgs, Nucl.Phys. B 261, 455 (1985). 

[23] T. Appelquist and H. D. Politzer, Phys. Rev. Lett. 34, 43 (1975). 

[24] C. Bernard, M.C. Ogilvie, T. DeGrand, C. DeTar, S. Gottlieb, A. Krasnitz, R. Sugar, 
and D. Toussaint, 'QCD thermodynamics with two flavors at Nt — 6' University of 
Arizona report AZPH-TH/91-45, (1991) 

[25] M.C. Chu, M. Lissia, and J.W. Negele, Nucl. Phys. B 360, 31 (1991). 

[26] M. Lissia et al, MIT CTP Preprint 1922 (1991). 

[27] G.E. Brown, M. Rho and W. Weise, Nucl. Phys. A 454, 669 (1986) . 

[28] F. Karsch, Z. Phys. C 38, 147 (1988). 

[29] D. Coffman et al., Phys. Rev. Lett. 68, 282 (1992). 



31 



[30] E. Shuryak, Nucl. Phys. B 214, 237 (1983). 

[31] E. Shuryak and J. L. Rosner, Phys. Lett. B 218, 72 (1989). 

[32] E. Shuryak, Nucl. Phys. A 533, 761 (1991) . 



32 



Figure captions: 



Figure 1: Schematic picture of replications leading to the two-dimensional Coulomb po- 
tential (|^) 



Figure 2: Wave function (full line) and potential (dashed line) for two-dimensional 
positronium with x = r{eT) 



Figure 3: Comparison of finite temperature Coulomb potentials. The solid lines cor- 
respond to the potential (|^), while the dashed ones correspond to the two-dimensional 
potential, 2Tlog(rT). 
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Figure 4: The potential of ref.pO[. The finite temperature (T > Tc) space-hke potential 
corresponds to the full line, the zero temperature time-like potential to the dashed line. 
The dashed-dotted line represents the finite temperature (T > Tc) space-hke potential and 
the short-dashed-long-dashed line the modified Coulomb potential according to eqn.(||). 



Figure 5: Wave function for the p-meson as function of r/a (a = 0.22 fm) for different 
temperatures. The upper set of curves corresponds to wave functions obtained with 
the modified Coulomb potential (eqn.([3^)) the lower set with the confining potential 
(eqn. (|3lD ). Temperatures: 150 MeV (full line), 250 MeV (short-dashed line) and 350 MeV 
(long-dashed line). Data are form ref. and correspond to T ^ 210 MeV. 



Figure 6: Wave functions for p and vr at T = 210 MeV (a) and T = 350 MeV (b) as a 
function of r/a {a = 0.23 fm). The full line and the short-dashed line correspond to the 
p with spin quantum number = ±1 and Sz = respectively. The long-dashed line 
represents the pion wave function. In (a) arbitrary normalization was used while in (b) 



all wave functions have the same normalization. Data are form ref. f 17 



Figure 7: Different paths contributing to a equal time correlator of two currents separated 
by a distance z. In the T = ) the transverse deviation of the paths is ~ 2; while 

at high T for non-interacting quarks (b) it is smaller, r± ~ \/z/T. The bound quarks (c) 

still deviate from the axis by r± ~ ^z/T, but their relative distance is defined by the size 
of the bound state rbound- 
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